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Abstract 

In invasion percolation, the edges of successively maximal weight 
(the outlets) divide the invasion cluster into a chain of ponds sepa- 
rated by outlets. On the regular tree, the ponds are shown to grow 
exponentially, with law of large numbers, central limit theorem and 
large deviation results. The tail asymptotics for a fixed pond are also 
studied and are shown to be related to the asymptotics of a critical 
percolation cluster, with a logarithmic correction. 

1 Introduction and definitions 

1.1 The model: invasion percolation, ponds 
and outlets 

Consider an infinite connected locally finite graph Q, with a distin- 
guished vertex o, the root. On each edge, place an independent 
Uniform[0, 1] edge weight, which we may assume (a.s.) to be all dis- 
tinct. Starting from the subgraph Co = {o}, inductively grow a se- 
quence of subgraphs Cj according to the following deterministic rule. 
At step i, examine the edges on the boundary of Cj_i, and form Cj 
by adjoining to Cj_i the edge whose weight is minimal. The infinite 
union 

oo 

C=\JCi (1.1) 

i=l 

is called the invasion cluster. 
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Invasion percolation is closely related to ordinary (Bernoulli) per- 
colation. For instance, ([3] for Q = Z d ; later greatly generalized by 
if Q is quasi-transitive, then for any p > p c , only a finite number 
of edges of weight greater than p are ever invaded. On the other hand, 
it is elementary to show that for any p < p c , infinitely many edges of 
weight greater than p must be invaded. In other words, writing £j for 
the weight of the z th invaded edge, we have 



So invasion percolation produces an infinite cluster using only slightly 
more than critical edges, even though there may be no infinite cluster 
at criticality. The fact that invasion percolation is linked to the critical 
value Pc, even though it contains no parameter in its definition, makes 
it an example of self- organized criticality. 

Under mild hypotheses (see section 13. 1|) , the invasion cluster has 
a natural decomposition into ponds and outlets. Let e\ £ C be the 
edge whose weight Qi is the largest ever invaded. For 72 > 1, e n is 
the edge in C whose weight Q n is the highest among edges invaded 
after e n -\. We call e n the n th outlet and Q n the corresponding outlet 
weight. Write V n for the step at which e n was invaded, with Vq = 0. 
The n th pond is the subgraph of edges invaded at steps i G (V n -x, V n ]. 

Suppose an edge e, with weight p, is first examined at step i G 
(y n -i,V n \. (That is, i is the first step at which e is on the boundary 
of Cj_i.) Then we have the following dichotomy: either 

• e will be invaded as part of the n th pond (if p < Q n ); or 

• e will never be invaded (if p ^> Qn) 

This implies that the ponds are connected subgraphs and touch each 
other only at the outlets. Moreover, the outlets are pivotal in the 
sense that any infinite non-intersecting path in C starting at o must 
pass through every outlet. Consequently C is decomposed as an infinite 
chain of ponds, connected at the outlets. 

In this paper we take Q to be a regular tree and analyze the 
asymptotic behaviour of the ponds, the outlets and the outlet weights. 
This problem can be approached in two directions: by considering the 
ponds as a sequence and studying the growth properties of that se- 
quence; or by considering a fixed pond and finding its asymptotics. We 
will see that the sequence of ponds grows exponentially, with exact ex- 
ponential constants. For a fixed pond, its asymptotics correspond to 
those of ordinary percolation with a logarithmic correction. 



limsup£i = p c 




oo 
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These computations are based on representing C in terms of the 
outlet weights Q n , as in [lj. Conditional on (Q n )^L , each pond is an 
independent percolation cluster with parameter related to Q n . In par- 
ticular, the fluctuations of the ponds are a combination of fluctuations 
in Q n and the additional randomness. 

Surprisingly, in all but the large deviation sense, the asymptotic 
behaviour for the ponds is controlled by the outlet weights alone: the 
remaining randomness after conditioning only on (Qn) ( ^ = Q disappears 
in the limit, and the fluctuations are attributable solely to fluctuations 
of Q n . 

1.2 Known results 

The terminology of ponds and outlets comes from the following de- 
scription (see |17j ) of invasion percolation. Consider a random land- 
scape where the edge weights represent the heights of channels between 
locations. Pour water into the landscape at o; then as more and more 
water is added, it will flow into neighbouring areas according to the 
invasion percolation mechanism. The water level at o, and throughout 
the first pond, will rise until it reaches the height of the first outlet. 
Once water flows over an outlet, however, it will flow into a new pond 
where the water will only ever rise to a lower height. Note that the 
water level in the n th pond is the height (edge weight) of the n th outlet. 

The edge weights may also be interpreted as energy barriers for a 
random walker exploring a random energy landscape: see [15]. If the 
energy levels are highly separated, then (with high probability and 
until some large time horizon) the walker will visit the ponds in order, 
spending a long time in each pond before crossing the next outlet. In 
this interpretation the growth rate of the ponds determines the effect 
of entropy on this analysis. See the extended discussion in [15]. 

Invasion percolation is also related to the incipient infinite cluster 
(IIC), at least in the cases Q = I? (|12j) and Q a regular tree: see, 
e.g., [32], [T] and [5]. For a cylinder event E, the law of the IIC can 
be defined by 

Piic(#) d = lim ¥ Pc (E I o «-► dB(k)) (1.3) 

k— >oo 

or by other limiting procedures, many of which can be proved to be 
equivalent to each other. Both the invasion cluster and the IIC consist 
of an infinite cluster that is "almost critical" , in view of (|1.2f) or (11.31) 
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respectively. For Q = I? ([12]) and Q a regular tree ([I]), the IIC 
can be denned in terms of the invasion cluster: if denotes a vertex 
chosen uniformly from among the invaded vertices within distance k 
of o, and T~x k E denotes the translation of E when o is sent to A% then 

P nc (£) = lim P{r Xk E) (1.4) 

fc^oo 

Surprisingly, despite this local equivalence, the invasion cluster and 
the IIC are globally different: they are mutually singular and, at least 
on the regular tree, have different scaling limits, although they have 
the same scaling exponents. 

The regular tree case, first considered in [16] , was studied in great 
detail in |1J. Any infinite non-intersecting path from o must pass 
through every outlet; on a tree, this implies that there is a backbone, 
the unique infinite non-intersecting path from o. In pQ a description of 
the invasion cluster was given in terms of the forward maximal weight 
process, the outlet weights indexed by height along the backbone (see 
section 13 . 2[) . This parametrization in terms of the external geometry 
of the tree allowed the calculation of natural geometric quantities, 
such as the number of invaded edges within a ball. In the following, 
we will see that when information about the heights is discarded, the 
process of edge weights takes an even simpler form. 

The detailed structural information in p] was used in [2] to identify 
the scaling limit of the invasion cluster (again for the regular tree). 
Since the invasion cluster is a tree with a single infinite end, it can be 
encoded by its Lukaciewicz path or its height and contour functions. 
Within each pond, the scaling limit of the Lukaciewicz path is com- 
puted, and the different ponds are stitched together to provide the full 
scaling limit. 

The two-dimensional case was also studied in a series of papers 
by van den Berg, Damron, Jarai, Sapozhnikov and Vagvolgyi (|17|. 
[6] and [5]). There they study, among other things, the probability 
that the n th pond extends a distance k from o, for n fixed. For n = 
1 this is asymptotically of the same order as the probability that a 
critical percolation cluster extends a distance k, and for n > 1 there 
is a correction factor (log/c)™ -1 . Furthermore an exponential growth 
bound for the ponds is given. This present work was motivated in 
part by the question of what the corresponding results would be for 
the tree. Quite remarkably, they are essentially the same, suggesting 
that a more general phenomenon may be involved. 
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In the results and proofs that follow, we shall see that the sequence 
of outlet weights plays a dominant role. Indeed, all of the results in 
Theorems 12.1 fj27il are proved first for Q n , then extended to other pond 
quantities using conditional tail estimates. Consequently, all of the re- 
sults can be understood as consequences of the growth mechanism for 
the sequence Q n . On the regular tree, we are able to give an exact de- 
scription of the sequence Q n in terms of a sum of independent random 
variables (see section I3T31) . In more general graphs, this representation 
cannot be expected to hold exactly. However, the similarities between 
the between the pond behaviours, even on graphs as different as the 
tree and Z 2 , suggest that an approximate analogue may hold. Such a 
result would provide a unifying explanation for both the exponential 
pond growth and the asymptotics of a fixed pond, even on potentially 
quite general graphs. 

1.3 Summary of notation 

We will primarily consider the case where Q is the forward regular tree 
of degree a: namely, the tree in which the root o has degree a and 
every other vertex has degree a + 1. The weight of the i th invaded 
edge is The n th outlet is e n and its edge weight is Q n . We may 
naturally consider e n to be an oriented edge e n = (v n ,v n ), where v n 
is invaded before v n . The step at which e n is invaded is denoted V n 
and the (graph) distance from o to v n is L n . Setting Vq = Lq = for 
convenience, we write V n = V n — V n -i and L n = L n — L n -\. 

There is a natural geometric interpretation of L n as the length 
of the part of the backbone in the n th pond, and V n as the volume 
(number of edges) of the n th pond. In particular V n is the volume of 
the union of the first n ponds. 

R n is the length of the longest upward-pointing path in the n th 
pond, and R' n is the length of the longest upward-pointing path in the 
union of the first n ponds. 

We shall later work with the quantity 5 n ; for its definition, see 
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We note the following elementary relations: 

n n 

L n = Y J U, V n = J2Vi, (1-5) 

i=i i=i 

n 

Qn+1 < Qn, < #n ^ (1-6) 

i=l 

Probability laws will generically be denoted P. For p £ [0,1], 
P p denotes the law of Bernoulli percolation with parameter p. For 
a set A of vertices, the event {x <-> ^4} means that there is a path 
of open edges joining x to some point of A, and {x <-> cxd} means 
that there is an infinite non-intersecting path of open edges starting 
at x. We define the percolation probability 9{p) = P p (o <-> oo) and 
p c = inf {p : 9{p) > 0}. dB(k) denotes the vertices at distance exactly 
k from o. 

For non-zero functions f(x) and g(x), we write /(#) ~ g(x) if 
lim = 1; the point at which the limit is to be taken will usually be 
clear from the context. We write f(x) x g(x) if there are constants c 
and C such that cg(x) < f(x) < Cg{x). 

2 Main results 

2.1 Exponential growth of the ponds 

Let Z n denote the 7-tuple 

Z n = ( log((Q n - PcY 1 ) , log L n , log L n , 

V . x (2-1) 

log R n , log R' n , \ log V n , \ log V n j 

and write 1 = (1, 1, 1, 1, 1, 1, 1). 
Theorem 2.1. PFii/i probability 1, 

lim — = 1. (2.2) 
Theorem 2.2. If (B t )t>o denotes a standard Brownian motion then 

( Z\ Nf \ — Nt ■ 1 \ 

-^7= (Bt ■ l)t>o (2-3) 

/ t>o 



N 



6 



as N — > oo ; with respect to the metric of uniform convergence on 
compact intervals oft. 

These theorems say that each component of Z satisfies a law of 
large numbers and functional central limit theorem, with the same 
limiting Brownian motion for each component. 

Theorem 12.21 shows that the logarithmic scaling in Theorem 12.11 
cannot be replaced by a linear rescaling such as e n (Q n — p c ). Indeed, 
log((<5n — Pc) -1 ) has characteristic additive fluctuations of order ±y/n, 
and therefore Q n — p c fluctuates by a multiplicative factor of the form 
e ±v/ ™. As n — > oo this will be concentrated at and oo, causing 
tightness to fail. 

Theorem 2.3. - log((Q n — Pc) 1 ) satisfies a large deviation principle 
on [0, oo) with rate n and rate function 

ip(u) = u — logn — 1. (2.4) 

MogL n , Mogi? n and ^logl^i satisfy large deviation principles on 
[0, oo) with rate n and rate function ip, where 



a. — logti — 1 ifu > 5, 



[log(2)-n »/«<5- 

It will be shown that ip arises as the solution of the variational 
problem 

tp(u) = inf (93(f) + v — u) (2-6) 

v>u 

2.2 Tail behaviour of a pond 

Theorems [27TH2T3] describe the growth of the ponds as a sequence. We 
now consider a fixed pond and study its tail behaviour. 

Theorem 2.4. For n fixed and e — * + , k — * 00, 

2a eHoge- 1 )"- 1 
P {Qn < Pc(l + e)) ~ — (2.7) 



and 



k) 
~~k 

¥(V n > k) x P [V n > k) x (1 ° g ^ W " 1 (2.10) 
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Using the well-known asymptotics 



2<7 2 

0(p)~ t(p-Pc) as p^p+, (2.11) 

a — 1 



2a 



\ c (o <-> dB(k)) ~ 7- — — as fc -> oo, (2.12) 



we may rewrite fl2.7p - fl2.10jl as 



1 (Q„ < p c (l + e)) ~ ^y.i), g (Pc(l + e)) (2.13) 



'(L„ >fc)~P(L n > A;) ~ \° ^, P Pe (o~aB(fc)) (2.14) 



(log A:)™" 1 
(n- 1)1 

P(R n > k) xP(i?; > k) x (logjb) n-1 Pp c (o^aB(Jfc)) (2.15) 
P(F n > k) xP (K > fe) x (log A;) n_1 P Pc ( | C(o) | > fc) (2.16) 

Working in the case Q = Z 2 , [6] considers i? n , the maximum dis- 
tance from o to a point in the first n ponds, which is essentially R' n in 
our notation. [6] Theorem 1.5] states that 

P(i?n > fc) x (log k) n ~% c (o «-► 8B{k)) (2.17) 

and notes as a corollary 

P(^>fc)xP Pc (o^3B(fc)) (2.18) 

where A denotes a percolation connection where up to i edges are 
allowed to be vacant ("percolation with defects"). (|2.18p suggests the 
somewhat plausible heuristic of approximating the union of the first 
n ponds by the set of vertices reachable by critical percolation with 
at most n — 1 defects. Indeed, the proof of (j2.17p uses in part a 
comparison to percolation with defects. By contrast, on the tree the 
following result holds: 

Theorem 2.5. For fixed n and k — ► oo, 

p Pc ( AaB(fc))xr r " (2.i9) 

The dramatic contrast between fl2.18H and (|2.19p can be explained 
in terms of the number of large clusters in a box. In Z 2 , a box of 
side length S has generically only one cluster of diameter of order S. 
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On the tree, by contrast, there are many large clusters. Indeed, a 
cluster of size N has of order N edges on its outer boundary, any one 
of which might be adjacent to another large cluster, independently of 
every other edge. Percolation with defects allows the best boundary 
edge to be chosen, whereas invasion percolation is unlikely to invade 
the optimal edge. 

2.3 Outline of the paper 

Section 13.11 states a Markov property for the outlet weights that is 
valid for any graph. From section 13.21 onwards, we specialize to the 
case where Q is a regular tree. In section 13.21 we recall results from 
PQ that describe the structure of the invasion cluster conditional on 
the outlet weights Q n . Section 13.31 analyzes the Markov transition 
mechanism of section 13.11 and proves the results of Theorems I2.1H2.3I 
for Q n . 

Section [4.11 states conditional tail bounds for L n , R n and V n given 
Q n . The rest of sections HHS use these tail bounds to prove Theorems 
I2.1H2.4I The proof of the bounds in section 14.11 is given in section [7J 
Finally, section [8] gives the proof of Theorem 12.51 

3 Markov structure of invasion perco- 
lation 

In section 13.11 we give sufficient conditions for the existence of ponds 
and outlets, and state a Markov property for the ponds, outlets and 
outlet weights. Section [3.21 summarizes some previous results from pQ 
concerning the structure of the invasion cluster. Finally in section 13.31 
we analyze the resulting Markov chain in the special case where Q is 
a regular tree and prove the results of Theorems 12.1142731 for Q n . 

3.1 General graphs: ponds, outlets and outlet 
weights 

The representation of an invasion cluster in terms of ponds and outlets 
is guaranteed to be valid under the following two assumptions: 

0(p c ) = (3.1) 
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and 

limsup£j=p c a.s. (3-2) 



(|3.ip is known to hold for many graphs and is conjectured to hold for 
any transitive graph for which p c < 1 ([3l Conjecture 4]; see also, for 
instance, [131 section 8.3]). If the graph Q is quasi-transitive, (|3.2p 
follows from the general result \\.\\ Proposition 3.1]. Both (|3.ip and 
(|3.2p hold when Q is a regular tree. 

The assumption (|3.ip implies that w.p. 1, 

sup£j > p c for all iq (3.3) 

i>io 

since otherwise there would exist somewhere an infinite percolation 
cluster at level p c . We can then make the inductive definition 

Qi = max& = &> (3.4) 
Q„ = max & = £^ (n > 1) (3.5) 



since (|3.2p and (|3.3|) imply that the maxima are attained. 

Condition on Q n , e n , and the union C n of the first n ponds. We 
may naturally consider e n to be an oriented edge e n = (v n ,v n ) where 
the vertex v n was invaded before v n . The condition that e n is an 
outlet, with weight Q n , implies that there must exist an infinite path 
of edges with weights at most Q n , starting from v n and remaining in 
Q\C n . However, the law of the edge weights in Q\C n is not otherwise 
affected by Q n ,e n ,C n . In particular we have 

n / 'In ^ ^\ P 9 '(e n ^ oo in Q\C n ) 

Qn+i < q \Qn,e n ,Cn) = ~ — . —< = : (3-6) 

' P(?„(en <-> oo m Q\C n ) 

on the event {q' < Q n }- In (|3.6p we can replace G\C n by the connected 
component of Q\C n that contains e n . 

3.2 Geometric structure of the invasion clus- 
ter: the regular tree case 

In [H section 3.1], the same outlet weights are studied, parametrized 
by height rather than by pond. is defined to be the maximum 
invaded edge weight above the vertex at height k along the backbone. 
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A key point in the analysis in pQ is the observation that (Wfc)S_ 
is itself a Markov process. Wk is constant for long stretches, corre- 
sponding to k in the same pond, and the jumps of Wk occur when an 
outlet is encountered. The relation between the two processes is given 
by 

W k = Q n iff L n _! <k<L n (3.7) 

From (|3.7p we see that the (Q n )^Lo are the successive distinct 
values of (Wfc)^L , and L n = L n — L n _i is the length of time the 
Markov chain Wk spends in state Q n before jumping to state Q n +i- 
In particular, L n is geometric conditional on Q n , with some parameter 
depending only on Q n . As we will refer to it often, we define 5 n to be 
that geometric parameter: 

F(L n >m\Q n ) = (l-5 n ) m (3.8) 

A further analysis (see [H section 2.1]) shows that the off-backbone 
part of the n th pond is a sub-critical Bernoulli percolation cluster 
with a parameter depending on Q n , independently in each pond. We 
summarize these results in the following theorem. 

Theorem 3.1 QdJ, sections 2.1 and 3.1). Conditional on (Qn)n°=i! 
the n th pond of the invasion cluster consists of 

1. L n edges from the infinite backbone, where L n is geometric with 
parameter 5 n ; and 

2. emerging along the a — 1 sibling edges of each backbone edge, 
independent sub-critical Bernoulli percolation clusters with pa- 
rameter 

Pc(l - S n ) (3.9) 

Given (Q n )^=o> the ponds are conditionally independent for different 
n. 5 n is a continuous, strictly increasing functions of Q n and satisfies 

5n ~ ^-Wn) ~ <y(Qn ~ Pc) (3.10) 

la 

as Q n -> pt- 

The meaning of (13.101) is that 5 n = f(Q n ) where f(q) ~ ^^"^(?) ~ 
cr(q-p c ) as q -> p+. 

It is not at first apparent that the geometric parameter 6 n in (13. 8|) 
is the same quantity that appears in (|3.9p . and indeed [Ij has two 
different notations for the two quantities: see [U equations (3.1) and 
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(2.14)]. Combining equations (2.3), (2.5), (2.14) and (3.1) of [I] shows 
that they are equivalent to 

5 n = l-<jQ n (l-Q n 6{Q n )y- 1 (3.11) 

For a = 2 we can find explicit formulas for these parameters: p c = i, 
Q(p) = p ~ 2 (2p - 1) for p > p c , 5 n = 2Q n - 1 and p c {\ - S n ) = 1 - Q n . 
However, all the information needed for our purposes is contained in 
the asymptotic relation (13,10p . 

3.3 The outlet weight process 

The representation (13. 6p simplifies dramatically when Q is a regular 
tree. Then the connected component of Q\C n containing e n is isomor- 
phic to G, with e n corresponding to the root. Therefore the depen- 
dence of Qn+i on e n and C n is eliminated and we have the following 
result. 

Corollary 3.2. On the regular tree, 
Markov chain with 

HQi < v) 

and 

P(Qn+l <q'\Qn 

for p c < q' < q. 

Equations (|3.12p and f)3. 13|) say that, conditional on Q n , Q n +i is 
chosen from the same distribution, conditioned to be smaller. In terms 
of (W / fc)^ =0 , (|3.13p describes the jumps of when they occur, and 
indeed the transition mechanism (|3.13p is implicit in pp. 

Since 8 is a continuous function, it is simpler to consider 9{Q n )'- 
0(Q\) is uniform on [0, 1] and 

7/' 

¥(e(Q n+1 )< U '\9(Q n )=u) = - (3.14) 

u 

for < u' < u. But this is equivalent to multiplying 8(Q n ) by an 
independent Uniform[0, 1] variable. Noting further that the logarithm 
of a Uniform [0, 1] variable is exponential of mean 1, we have proved 
the following proposition. 



(Qnj'^Li is a time-homogeneous 
= 9{q) (3.12) 

-4-m {3A3) 
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Proposition 3.3. Let Ui, i S N ; be independent Uniform[0, 1] random 
variables. Then, as processes, 




(3.15) 



Equivalently, with Ei = log U i independent exponentials of mean 1, 

n 

log^Qn)" 1 ) =J>i (3.16) 

i=l 

jointly for all n. 

Corollary 3.4. The triple 

Z' n = (log(0(Q n )- 1 ),log((g n - Pc )- 1 ),log ( 5- 1 ) (3.17) 

satisfies the conclusions of Theorems \2.1\ and \2.2\ and each compo- 
nent of \Z' n satisfies a large deviation principle with rate n and rate 
function 

ip(u) = u — logu — 1. (3.18) 

Proof. The conclusions about log (9(Q n )^ 1 ) follow from the repre- 
sentation (|3.16p in terms of a sum of independent variables; the rate 
function ip is given by Cramer's theorem. The other results then follow 
from the asymptotic relation (|3. 10[) . □ 

4 Law of large numbers and central 
limit theorem 

4.1 Tail bounds for pond statistics 

Theorem 13.11 expressed L n ,R n and V n as random variables whose pa- 
rameters are given in terms of Q n . Their fluctuations are therefore 
a combination of fluctuations arising from Q n , and additional ran- 
domness. The following proposition gives bounds on the additional 
randomness. 

Recall that S n is a certain function of Q n with 5 n ~ o~(Q n — p c ): 
see Theorem 13.11 
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Proposition 4.1. There exist positive constants C, c, so, 7l, 7_r, "fv 
such that L n , R n and V n satisfy the conditional bounds 

V{5 n L n > S\5) <Ce~ cS P(S n L n <s\5) <Cs (4.1) 
W(6 n Rn > S\5) < Ce~ cS V(5 n Rn< s\S) <Cs (4.2) 
P(^y n > S\ 5) < Ce~ cS F(5lV n < s\S) < C^fs (4.3) 

for all n and all S,s > 0; and 

P ( 5 n L n <s\S)>cs on {5 n < 7L s} (4.4) 

P ( SnRn <s\5)>cs on {5 n < lR s} (4.5) 

W(5lV n <s\6) >c^~s on{5 2 n < lv s} (4.6) 

for s < s - 

The proofs of (|4.ip ~ (|4.6p . which involve random walk and branch- 
ing process estimates, are deferred to section [71 

4.2 A uniform convergence lemma 

Because Theorem 12.21 involves weak convergence of several processes to 
the same joint limit, it will be convenient to use Skorohod's represen- 
tation theorem and almost sure convergence. The following uniform 
convergence result will be used to extend convergence from one set of 
coupled random variables (<5 n ,7v) to another (X n ^): see section [4T3l 



Lemma 4.2. Suppose {X Ut N} n Ne fq and {5 n ^} n NeN are positive ran- 
dom variables such that 5 n ^ is decreasing in n for each fixed N, and 
for positive constants a, j3 and C , 

n8 a n>N X n>N >S)< CS-? (4.7) 
P(8* >N X n , N <s)< Cs? (4.8) 

for all S and s. Define 

n 

X ntN = Y J X l , N . (4.9) 
i=l 

Then for any T > and a > 0, w.p. 1, 

\og{5l N X n , N ) \og{5l N X n , N ) 
iim max = lim max = (J. 

TV^oo l<n<NT N a N^oo l<n<NT N a 

(4.10) 
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Proof. Let e > be given. For a fixed N, (|4.7p implies 

lo g(Ci\r^n,iv) 



max 



Kn<NT N a / ~ ^ V n,W ' 

_ _ / l<n<NT 

n ( e Nae \ 

< E E P (<NXi,N>—) 

l<n<NT i=l V ' 

n f p JV°e\ 

< E E p w>v 



Kn<NT i=l 



< E E^" 



-/JAT a e 
l<n<NT i=l 

< (NTf^Ce-^ (4.11) 



where we used > 5 nt N in the third inequality. But then since 
Y^N=i(NT) 2+ PCe~P Na>L < oo, the Borel-Cantelli lemma implies 

1- Mff, A*) 

limsup max — < e 4.12 

a.s. Similarly, (|4.8p implies 

Pf max 1 ° g( ^ X "' Af) < < V P(^ jV X nJV <e-^) 
v - - 7 \<n<NT 

< NTCe- pNae (4.13) 

so that 

limmt max > — e (4.14) 

iV^oo \<n<NT N a 

a.s. Since e > was arbitrary and X n ^ < X,iv, (|4.1U|) follows. □ 



4.3 Proof of Theorems EZEH2J2 



The conclusions about Q n are contained in Corollary 13.41 The other 
conclusions will follow from Lemma 14.21 From Corollary 13.41 we may 
apply Skorohod's representation theorem to produce realizations of 
the ponds for each N £~N, coupled so that 



N 



I (B t )o<t<T (4.15) 

0<t<T 
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cl.S. ELS N — > oo. Then the relation 



± log X niN -Nt _ lQ g (^L^J ,n) ~ Nt lQ g ( 5 f^J ,N X lNt\ ,n) 

(4.16) 

shows that = log X n will satisfy a central limit theorem as well. The 
same holds for X n . We will successively set 



X n ,N 


= L n> N, 


with a = 


1, 


(4.17) 


X ny N 


= Rn,N, 


with a = 


1, 


(4.18) 


X n ,N 


= V nt N, 


with a = 


2. 


(4.19) 



The bounds (|4.7|) - (|4.8|) follow immediately from the bounds in Propo- 
sition 2J3 This proves Theorem 12 .21 for L n and V n . For R n , the quan- 
tity R is not the one that appears in Theorem 12.21 but the bound 
Rn < R' n — Rn implies the result for R' n as well. 

The lemma also implies the law of large numbers results (I2.2p . by 
taking T = 1 and using the same ponds for every N. 

5 Large deviations: proof of Theorem 

EH 



In this section we present a proof of the large deviation results in 
Theorem l2.31 As in section^ we prove a generic result using a variable 
X n and tail estimates. Theorem 12.31 then follows immediately using 
Corollary 13.41 and Proposition 14.11 

Note that Proposition 15.11 uses the full strength of the bounds in 
Proposition 14.11 

Proposition 5.1. Suppose that 5 n and X n are positive random vari- 
ables such that, for positive constants a, /3, c, C, 7, sq, 

V(%X n > S\8 n )<Ce- cSf> (5.1) 

¥(8 a n X n <s\5 n )<Cs l / a (5.2) 

for all S and s, and 

¥{5 a n X n <s\5 n )>cs 1 l a (5.3) 
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on the event {<5° < 7s}, for s < sq. Suppose also that ^ logS^ 1 satis- 
fies a large deviation principle with rate n on [0, 00) with rate function 
ip such that 99(1) = 0, if is continuous on (0, 00), and ip is decreasing 
on (0,1] and increasing on [l,oo). Then -^logX n satisfies a large 
deviation principle with rate n on [0, 00) with rate function 

ip(u) = inf (ip(v) + v — u) (5-4) 

v>u 

Proof. It is easy to check that ip is continuous, decreasing on [0, 1] 
and increasing on [l,oo), ^(1) = 0, and ip{u) = <p{u) for u > 1. So it 
suffices to show that 

1/1 \ 
lim — logP — logX„ > it = — inf ip(v) (5.5) 
n— too n \ an I v>u 



for u > and 



lim - logP ( — logX n < u ) = -Mu) (5.6) 
n-*oo n \an 



for < u < 1. For (jo3|) . let e > 0. Then 

— logX„>u) <P[-logS- 1 >u-e 
an I \ n 



+ P ( ilog^ 1 < u -e, — logX„ > u 
\n an 

< P (- logd- 1 > u - e ) + P ( — log(%X n ) > e 
\n J \an 

< P log > u - + Ce~ cePant (5.7) 

The last term in (|5.7p is super- 
exponentially small, so (|5.7p and the large deviation principle for 
ifog^ 1 imply 

lim sup — log P ( — log X n > u J <— inf tp(v). (5.8) 
n— »oo n \an / v>u-e 

On the other hand, 

— log X n > u J > P f - log > « + e, — log(^X n ) > -e ) 
an J \n an J 

^pf-logS" 1 >u + ej (l-Ce~ ne ) (5.9) 
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using (|5.2p with s = e ane . So 

liminf — logP ( — logX n > u ) > - inf (f(v). (5.10) 



n^oo n \ an I v>u+e 



Since tp is continuous and e > was arbitrary, this proves (|5.5p . 

For (|5.6p . let u € (0, 1) be given and choose v £ (u, 1), e G (0, it). 
Then for n sufficiently large we have 

P (— logX n < u) > P (v - e < - log J" 1 < v, — log(^X n ) <u-v) 
\an J \ n an J 

> P (v - e < i log (J" 1 < ^ ce-™^-") (5.11) 

Here we used (I5.3P with s = e _an ( v_u ) . Note that if n is large enough 
then s < so and the condition 5® < js follows from v — e < - \og8~ 1 . 
Therefore, since ip is decreasing on (0,1], 

liminf — logP [ — logX„ < u J > — ( inf <p{w) J — (u — it) 



n— >oo n \an / \y~e<w<v 

= -tp(v) - a(v - u) (5.12) 

(|5.12p was proved for u < v < 1. However, since </? is continuous and 
the function —<p(v) — av is decreasing in v for u > 1, (|5.12p holds for 
all v > u. So take the supremum over v > u to obtain 

liminf -logP f — log X n < u J > -^(u). (5.13) 

Finally 
1 



log X n < w 

an 



< 



-log^ 1 < u) +E (1 f -log^ 1 > u J P ( — logX n < u 
n / \ \n J \an 



- log 5 n 1 < u 

n 



+ pflf-log(5- 1 >^P( — log(<5£X n ) <u --log5- 1 
^log^ 1 < u \ +E (l (^log^n 1 > A Ce"^-^ 1085 " 1 )^ 



< _ 

n 



(5.14) 
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(using (pT2j) with s = e an ( u ™ ^s^ 1 )). Apply Varadhan's lemma (see, 
e.g., p. 32]) to the second term of f)5 . 14|) : 

lim ilogE (l ( -\ogS~ 1 > u) Ce n ( u -^ losS ^)) 
n^oon \ \n J J ( 515 ) 

= sup(u — v — f(v)) = —ip(u) 

v>u 

Therefore 

1/1 \ 
lim sup — logP — logX n < u < max{— <p(u), —tp(u)} = —tp(u) 
n^oo n \an J 

(5.16) 

which completes the proof. □ 



6 Tail asymptotics 

In this section we prove the fixed-pond asymptotics from Theorem [27 



Proof of d22|). Recall from fl3~T6l) that \og(9{Q n )~ l ) has the same 
distribution as a sum of n Exponential variables of mean 1, i.e., a 
Gamma variable with parameters n, 1. So 

P(0(Qn) < e ) = P (log(^Qn)" 1 ) > lege" 1 ) 

Make the substitution x = (1 + u) loge -1 : 

/i —IN™ r-cx) 

P(0(Q n ) <e) = € \° S€ J / (1 + uT^e-^' 1 (6.2) 



(n-1)! .„, 

Then Watson's lemma (see for instance (2.13) of p3|) implies that 

>n-l 



n0{Qn)<e)^ e{lOge ~ 1 } ((U) 
(n — 1)! 



and so (|3.10p gives 



pw„<p e (i + .))~ 2 " (i °f (6.4) 

(a- l)(n- 1)! 

□ 
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Combining ()6.3[) with (|3. 10[) implies at once that 

P(«5 n <e)xe(loge- 1 )"- 1 (6.5) 
We use f)6.5[) to prove (|2.9p - (|2.10p using the following lemma. 

Lemma 6.1. Let 5 n be a random variable satisfying (|6.5p . Suppose 
a, {3 are positive constants such that a/3 > 1, and X n is any positive 
random variable satisfying 

P(5%X n > S\5) <CS~ P (6.6) 

for all S,n > 0, and 

F(S%X n >s \S)> Po (6.7) 
for some so,Po > 0. Write X n = Y^i=\ X%. Then 

F(X n >^P(1„>^ {l ° g k ^ 1 (6.8) 

as k — > oo. 

Proof. From (|6.7p and (|6.5p . 

P(X n > fc) > p (a»x n > s 1 5 n < s fc~ 1/o ) P(5 n < so^" 1/a ) 
(logfc)"- 1 

For the lower bound we condition on the Gamma random variable 
log0(Q n ) -1 . From (|3.1Up we have 5 n > C20(Q n ) for some constant 
C2 > 0, so that 

F(X n >k)=E (P (8 a n X n > H a n | b n )) 

<e(i a c(k5«yP 

< C 3 E (l A {kQ(Q n ) a )~ P 
= C 4 y ^1 A (£;e~ aa: ) J x n ~ 1 e~ x dx (6.10) 
Make the substitution y = ke~ ax to obtain 



[X n >k)<C 5 J^ (l A y-P) (log fc - log yY 



C b {\ogk) n - 1 f°°, / , |logy|\ n_1 



< ,°/ / 1AjT p U + T^f -FT/" ( 6 - n ) 

fci/a y n V / V log A; J yl-l/a v 7 



20 



The last integral in (|6.1ip is bounded as k — > oo: the singularity as 
y —* + is integrable since 1 — 1/a > 0, and the singularity as y — > oo 
is integrable since the exponent in y- 1 ~P+ 1 / a has — 1 — (3 + l/a < — 1 
since a/3 > 1. 

To extend (|6.8p to X n , assume inductively that P(X n > k) x 
(logfc) n_1 A 1/a - (The case n = 1 is already proved since Xi = X\.) 
The bound P(X n+ i > k) > P(X n+ i > fc) is immediate, and we can 
estimate 

> fc) < P (X n+ i >k-k')+F (X n > k') (6.12) 



where we set k' = [A; / (log A;) a / 2 J . Then k — k' ~ k and log (A; — k') ~ 
log k' ~ log k, so that 

P (X n+1 > k - k') x (6.13) 

while 

which is of lower order. This completes the induction. □ 



Proof of (|2.9p - (j2.10p . These relations follow immediately from (|6.5p 
and Lemma I6.lt the bounds (I6.6p - (l6.7p are immediate consequences 
of Proposition 14.11 As in section 14.31 the asymptotics for R' n follow 
from the asymptotics for R n and the bound R n < R' n < R n . □ 

Proof of ([2SD. For L n , we can use the exact formula P (L n > k | (5^) — 
(1 — 5 n ) k . Write 5 n = g(0(Q n )), where g(p) is a certain continuous 
and increasing function. By (|3,10p . g(p) ~ ^^r-P as p — > + ; as above 
we will use the bound g(p) > cp for some constant c > 0. Proceeding 
as in (|6.10p and (|6.1ip we obtain the exact formula 

P(Ln >k) = fc(n-l)! J^ 1 ~ 9{y/k)) k {\ogk- logy)"" 1 dy 



k(n - 1)! J 



(6.15) 



But the integral in (|6.15p converges to J °° e~^~ y dy = as k — > oo: 
pointwise convergence follows from g{p) ~ ^^-P> an d we can uniformly 
bound the integrand using 

(1 - g(y/k)) k < e -W=) < e~ cy (6.16) 
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Lastly, a simple modification of the argument for X n extends (|2.8p to 
L n . □ 



7 Pond bounds: proof of Proposition 

am 

In this section we prove the tail bounds (14,lj) - (|4.6p . Since the laws of 
L n , R n and V n do not depend on n except through the value of S n , 
we will omit the subscript in this section. For convenient reference we 
recall the structure of the bounds: 

¥(5 a X > S\S) < Ce~ cS , (7.1) 

F(5 a X < s\S) < Cs (7.2) 

for all S and s, and 

P ( 5 a X < s\6) > cs 1/a (7.3) 

on the event {5 a < 7 s}, for s < sq. We have a = 1 for X = L and 
X = R, and a = 2 for X = V. 

In (|7.3p it is necessary to assume 5 a < 7s. This is due only to 
a discretization effect: for any N- valued random variable X, we have 
P (5 a X < s I 5) = whenever 5 a < s. Indeed, the bounds (|44l)- (f46l) 
can be proved with 7 = 1, although it is not necessary for our purposes. 

Note that, by proper choice of C and sq, we can assume that S 
is large in (I7.ip and s is small in (]7.2p and (|7.3I) . Since we only con- 
sider N- valued random variables X, we can assume 8 is small in (|7.2p . 
say 5 < 1/2 (otherwise take s < (l/2) a without loss of generality). 
Moreover, Theorem 13.11 shows that L, R and V are all stochastically 
decreasing in 5. Consequently it suffices to prove (17. ip for 5 small, say 
5 < 1/2. Finally the constraint 5 a < jsq makes 5 small in (|7.3p also. 

We note for subsequent use the inequalities 

(1 - x) y < e~ xy (7.4) 

for x S (0, 1) , y > and 

1 - (1 - x) y < 1 - e~ 2xy < 2xy (7.5) 

for x G (0, 1/2), y > 0, which follow from log(l — x) < —x for x G (0, 1) 
and log(l -x)> -2x for x G (0, 1/2). 



22 



7.1 The backbone length L: proof of ( 14.11 ) and 

(S3D 

From Theorem 13.11 L is a geometric random variable with parameter 
5. So 

P(L > S/S | 5) = (l-5) [s/si 

< e sis/s\ < e -s+5 < e -s+i (7 6) 

since <5 < 1, proving (|7.1|) . For (|7.2|) and (|7.3|) . we have 

F(L<s/5\5) = l-(l-6)^^- 1 . (7.7) 

For 5 < I we can use (I7.5|) to get 

P(L < s/5| <5) < 2£([s/<5] - 1) < 2s (7.8) 

which proves (|7.2p . For (17. 3p . take 7^ = 1/2. Then on the event 
{5 < 7ls} we have \s/5~\ > 3 so that expanding (|7.7|) as a binomial 
series gives 

F(L<s/5\5)>( \s/5] - 1) 5 - \ ( fs^ - 1) ( f^] - 2) 5 2 
>(s-<5)-^(s-5) 

^-^i (7 - 9) 

for s < 1 = s - So <|73J holds. 

7.2 The pond radius i?: proof of (D) and (EH) 

Conditional on <5 and L, R is the maximum height of a percolation 
cluster with parameter p c (l — 5) started from a path of length L. 
We have R > L so (j7.2[) follows immediately from the corresponding 
bound for L. R is stochastically dominated by 

L + max Ra (7.10) 

l<i<L 

where Ri is the maximum height of a branching process with offspring 
distribution Binomial(cr,p c (l — 5)) started from a single vertex, inde- 
pendently for each i. Define 



a k = P (Ri > k 



(7.11) 
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for k > 0. Thus a k is the probability that the branching process 
survives to generation k + l. By comparison with a critical branching 
process, 

Cl 



a-k < 



k 



for some constant C\. On the other hand, a k satisfies 

a-k+i = 1 - /(I - a-k) 



(7.12) 



(7.13) 



where _f(z) is the generating function for the offspring distribution 
of the branching process. (This is a reformulation of the well-known 
recursion for the extinction probability.) In particular, since f'(z) < 
/'(l) = a Pc (l -5) = 1-5, 

a k+1 < f(l)a k = a k (l-5). (7.14) 

Combining (f7TT2|) with ([7TT4"1) . 

a k+j < a k (l - Sy < ^-e- 5 i (7.15) 
and taking k = \S/25] > S/25, j = [S/5\ - \S/25] > S/25 - 2, 



d 



< 



< 



a is/s\ - s 
C 2 6e- C3S 



2C\b s(S/28~2) 



S 



(7.16) 



Using this estimate we can compute 
S 



R > 



< 



r S 

L> Ts 



S ) + P I L < — , Ri > — for some i < S/28 
2d 2d 



< C 4 e- C55 
Similarly 



+ [ Y5 



S \ Ca5e- C ^ s 



S 



<C 8 e 



-cgS 



6 

(7.17) 



5) > 



L< i'^ < ^ f0ralH< i 



> c w s 1 



C u 5e- Cl2S ^ 3/25 



> ci 3 s 



(7.18) 



provided 5 is sufficiently small compared to s, i.e., provided 'Jr is small 
enough. 
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7.3 The pond volume V: proof of (14.31 ) and 

(USD 

From Theorem 13. 11 conditional on 8 and L, V is the number of edges 
in a percolation cluster with parameter p c (l — 8), started from a path 
of length L and with no edges emerging from the top of the path. We 
can express V in terms of the return time of a random walk as follows. 

Start with an edge configuration with L backbone edges marked 
as occupied. Mark as unexamined the (a — 1)L edges adjacent to the 
backbone, not including the edges emerging from the top. At each 
step, take an unexamined edge (if any remain) and either (1) with 
probability 1 — p c (l — 8), mark it as vacant; or (2) with probability 
p c (l — 8), mark it as occupied and mark its child edges as unexamined. 
Let Nk denote the number of unexamined edges after k steps. Then it 
is easy to see that is a random walk Nf. = No + Yli=i Yi ( a * least 
until Nk = 0) where Nq = (a — 1)L and 



Yi U-l ™*-f-^ (7.19) 
[-1 w.p. 1 -p c (l - 8). 

Let T = M{k : N k = 0}. (T is finite a.s. since E (Y { \ 6) = -S < 
and Nk can jump down only by 1.) T counts the total number of off- 
backbone edges examined, namely the number of non-backbone edges 
in the cluster and on its boundary, not including the edges from the 
top of the backbone. Consequently 

T = [V - L] + [(a - 1)V + a] - a = aV - L (7.20) 

In order to apply random walk estimates we write Xi = Yi + 5, = 
ELi X ^ so that E (Xi | 8) = 0; ci < E (Xf \ 5) < C 2 for universal 
constants c\, C2; and Nk = Zj~ — k5 + (a — 1)L. Note that 

{V > V } = {T> aV - L} 

C {Z laVo _ L} > 5 [aV -L\-(a- 1)L} (7.21) 

so using, for instance, Freedman's inequality p[ Proposition 2.1] leads 
after some computation to 

P ( V > S L/8 \8,L)<¥( Z [L{(TSo _ s y si > L(aS - 28 - a + 1) | 8, L) 

( L 2 {aS -28-a + l) 2 

- eXP V 2(aL(aS - 28 - a + 1) + C 2 L(aS - 8) /S) 
< exp (-c 3 S L8{l - 2/S ) 2 ) < exp(-c 4 S L8) (7.22) 
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if So > 3, say. Then, setting So = S/5L, 

P (V > S/5 2 | 5) < P (L > S/35 \5)+F(V > S L/5, L < S/35 \ 5) 

< C 5 e- C ° s + exp {-c A (S/5L)L5) 

< C 7 e~ CsS 

which proves (|7.ip . 

For (|7.2p . apply Freedman's inequality again: 

F(V < s/5 2 \5,L) =P(T <as/5 2 - L\S,L) 



(7.23) 



min (Zu — k5) < —(a — 1)L 



5,L 



<P{ min Z k < -(L- 5(as/5 2 )) 

\ k<crs/5 2 

( (L-as/5) 2 
" eXP V 2{a{L-as/5) + C 9 s/5 2 ) 
< exp (—c(SL — as) 2 / s) 



5,L 



(7.24) 
(7.25) 



(where in the denominator of (17.24p we use L < s/5 2 since V > L). 
So 



(V<s/5 2 \5) < F(L < 2as/5\S) 

+ P ( 1(l>2<j S /5) exp (-c(5L - as) 2 /s) \ 5) 



(7.26) 



The first term in (|7.26p is at most Cios from (|4.ip and will therefore 
be negligible compared to y/s. For the second term, note that for 
L > as/5, 



exp (— c(5L — as) 2 / s) 



crs/8 



l {l>L) 2cm ~ aS) e-^-^ 2 ^dl 



and so, with the substitution xy^ = 51 — as, 

E ( l(L>2as/5) exp {-c(5L - asf/s) \ 5) 



(7.27) 



2c 
2c 



' (L <l\5) 



6{5l - as) 



V <* 



{ be cx dx 



r 00 _ 2 

< Cn / {x^fs + crs)xe ca ' <ix < Ci2\/s 



(7.28) 
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which proves (|7.2[) 

Finally, for (|T.3[) . the Berry-Esseen inequality (see for instance [HJ 
Theorem 2.4.9]) implies 



[z k < -xVk^E (X 2 | 5) ^ - $ 



(-*) 



- t! (7,29) 



where <3?(x) = P(G < x) for G a standard Gaussian, and C13 is some 
absolute constant. In particular, using < c\ < E (X 2 | 5), 



(z k < -(a-l)Vk 



5) > C14 > 



(7.30) 



for k > C15. Choose = 1 A 7^ A (C15 + 1) . Then we have 
s/8 2 > 1 (so we may bound as/5 2 — yfs/5 > s/5 2 ); yfs/5 < 7l; and 
[s/5 2 \ > C15, so that 

P (V < s/5 2 I 5) = P (T < as/5 2 — L | 5) 

> P (T < as/5 2 - v^A L < y/s/6 \ 5) 

> P(T < s/5 2 ,L< x/s/5\6) 

> P ( Z Ls/52J < 5 [s/<5 2 J - (a - 1)L, L < ^fs/5 I 5) 

> p(z Ls/52 j < -(a- l) v r sJP^5^¥(L < yfi/6\6) 

> cie^ (7.31) 



proving (|7.3p . 



8 Percolation with defects 

In this section we prove 

F Pc (oAa#))xr r " (8.1) 

The case n = is a standard branching process result. For n > 
0, proceed by induction. Write C(o) for the percolation cluster of 
the root o. The lower bound follows from the following well-known 
estimate: 

P Pc (|C(o)|>iV) xi v-i/2 (8 .2) 

If C(o) > N then there are at least N vertices v\, . . . ,vjy on the outer 
boundary of C(o), any one of which may have a connection to dB(k) 
with n — 1 defects. As a worst-case estimate we may assume that 
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ui, . . . , ujv are still at distance k from dB(k), so that by independence 
we have 

F Pc (o A dB(k)) > P Pc (|C(o)| > iV) (l - (l - P Pc ( G 'S. 1 dB(k))^ ^ ^ 

>^(i-(i- C2fe -^ + y) (8.3) 



for constants ci , C2 . If we set N = k 2 " then the second factor is of 
order 1, and the lower bound is proved. For the upper bound, use a 
slightly stronger form of (|8.2j) (see for instance [lOl P- 260]): 

P Pc (|C( O )|=A0x(iV+ir 3 / 2 (8.4) 

Now if C(o) = N, with N < k/2, then there are at most o~N vertices 
on the outer boundary of C(o), one of which must have a connection 
with n — 1 defects of length at least k — N > k/2. So 

¥ Pc (o^dB(k)) 

<F Pc (\C(o)\>k/2) 

+ J2^ Pc (\C(o)\ = N)(l-{l-F Pc (0 n ^dB(k/2))y ) 

oo 



Ar=0 



< Ck^ 1 + c 7 (k 2 - n+1 Y 112 (8.5) 



which proves the result (the first term is an error term if n > 2 and 
combines with the second if n = 1). 
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